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Abstract 

We introduce a parametrically driven Ginzburg-Landau (GL) model, which 
admits a gradient representation, and is subcritical in the absence of the 
parametric drive (PD). In the case when PD acts uniformly in space, this 
model has a stable kink solution. A nontrivial situation takes places when 
PD is itself subject to a kink-like spatial modulation, so that it selects real 
and imaginary constant solutions at x — ±oo. In this situation, we find 
stationary solutions numerically, and also analytically for a particular case. 
They seem to be of two different types, viz., a pair of kinks in the real and 
imaginary components, or the same with an extra kink inserted into each 
component, but we show that both belong to a single continuous family of 
solutions. The family is parametrized by the coordinate of a point at which 
the extra kinks are inserted. However, solutions with more than one kink 
inserted into each component do not exist. Simulations show that the former 
solution is always stable, and the latter one is, in a certain sense, neutrally 
stable, as there is a special type of small perturbations that remain virtually 
constant in time, rather than decaying or growing (they eventually decay, 
but extremely slowly). 
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1 Introduction 



Various forms of parametrically driven Ginzburg-Landau (GL) equations con- 
stitute a class of pattern-forming systems with broken intrinsic phase-rotation 
invariance Jl|, @ . In this work, we consider a gradient model of this type, 

SL 

u t = e(x) ■ u* + (u + u xx - \u\ 2 u = -— , (1) 

ou* 

where u* stands for the complex conjugation, S/Su*is the variational deriva- 
tive, the function e{x) and constant ( are real, and the Lyapunov functional 
is 

L = - [-e(x) • (u 2 + (u*) 2 ) - (\u\ 2 + \u x \ 2 + \u\ 4 ]dx. (2) 
2 7— oo ^ ' 

Equation (1) with e = const has various applications, the best-known one be- 
ing the Rayleigh-Benard convection in the case when the convection-driving 
temperature difference across the fluid layer is subject to a resonant periodic 
modulation along the layer (the coordinate x is running along the layer), 
or the boundary is similarly undulated, see early || and more recent [|J 
theoretical works on this topic, and an experimental work ||. The "reso- 
nant" character of the modulation means that its period is half the depth 
of the convective layer. The presence of the variable coefficient e(x) implies 
that the "primary" modulation is further subject to a "supermodulation" on 
a scale essentially larger than the layer's depth, which is quite possible in 
large-aspect-ratio convection cells, and gives rise to novel patterns, as it will 
be shown below. 

The real constant ( in Eq. (1) controls the system's overcriticality in the 
absence of the parametric drive: when e = 0, the situation is oweroritical 
if ( > 0, and stzfrcritical if ( < 0. Rescaling allows one to reduce ( to one 
of the three values, ( = —1,0, or + 1. To the best of our knowledge, the 
previous works were dealing with the case when the gradient GL system 
was overcritical when e = 0. Our aim is to consider the opposite case with 
C = — 1, when no nontrivial pattern may exist without the parametric drive. 
Thus, C = — 1 will be assumed below, unless another value of ( is specified. 

In the case e = const, Eq. (1) has obvious constant solutions, 

Mo = ±V e — 1, or u = ±W— e — 1, (3) 
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which exist provided that, respectively, e > 1 or e < -1, while the trivial 
solution u = is stable in the case |e| < 1 , when the nonzero solutions (H) 
do not exist. In the case |e| > 1, another available exact solution is a kink, or 
domain wall, which provides for a transition between the constant solutions 
(HD with the opposite signs. For instance, in the case e > 1 (recall we set 
C = —1), the kink solution to Eq. (1) is 



u k mk{x) = ay/e- ltanh ^y — I > ( 4 ) 

where a = ±1 is the kink's polarity. 

Note that a similar kink, with (e — 1) replaced by 1 is, simultaneously, an 
exact solution to the usual GL equation, which can be obtained from Eq. (1) 
setting e = and ( = +1. It is commonly known that, in its latter capacity, 
the real kink solution is unstable against purely imaginary perturbations. 
On the contrary to this, an elementary consideration demonstrates that the 
solution (U) is stable as a solution to the parametrically driven GL equation 
(1) with ( = — 1 and e > 1. This is a formally new result, but the actual 
objective of the work is to consider more interesting solutions in the case 
when the parametric drive is ramped, i.e., it is itself subject to a kink-like 
modulation, see below. 



2 The ramp and stationary solutions 

A configuration with e(x) monotonically varying between asymptotically con- 
stant values at x = ±oo is called a ramp. In the case of the usual (not para- 
metrically driven) GL equation, much work has been done for a case when 
the overcriticality is ramped so that the system is subcritical at x = — oo and 
supercritical at x — +oo, as this configuration makes it possible to uniquely 
solve the wavenumber- selection problem, and demonstrates other remarkable 
features 0. Moreover, the ramp problem can be extended into the two- 
dimensional geometry, replacing u xx by the two-dimensional Laplacian. In 
particular, the stability of the ramp-supported quasi-one-dimensional pattern 
in the two-dimensional version of the usual GL equations was investigated in 
Refs. ||. In the parametrically driven equation (1) with £ = — 1, a similar 
ramp can easily be arranged too, choosing a function e(x) which is monoton- 
ically varying between a value |e(— oo)| < 1, that provides for the stability of 
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the u = solution at x = — oo, and |e(+oo)| > 1, which gives rise to a stable 
nonzero solution @ at x — +00. 

The subject of this work is to consider a more nontrivial case of the 
ramp, when e(x) interpolates between values e(— 00) < —1 and e(oo) > +1, 
for instance 

e(x) = a tanh (Ax) , a > 1, A > 0. (5) 

This ramp supports purely imaginary and purely real nonzero states (|J) 
at x = —00 and x = +00, so that u(x = —00) = ±iy/a — 1 , u(x = 
+00) = ±V fl — 1- Our objective is to find solutions interpolating between 
these asymptotic states, and investigate their stability. A configuration of 
the present type can be easily implemented experimentally, which makes 
the problem physically relevant (in the convective layer, the constant imagi- 
nary and real solutions represent uniform arrays of convection rolls mutually 
shifted by half the layer's depth). 

For the numerical solution of Eq. (1) with ( = —1 and e(x) taken as per 
Eq. (P), we define u(x) = v(x) + iw(x), arriving at a system of equations for 
two real functions, 

v t = a tanh (Ax) • v + v xx — (v 2 + w 2 + lj v = , (6) 
w t = —a tanh (Ax) • w + w xx — (v 2 + w 2 + w = . (7) 

Boundary conditions (b.c.) corresponding to the situation defined above are 



v(x = —00) = 0, v(x = +00) = ±ya — 1, (8) 

w(x = —00) = ±\/a — 1 , w(x = +00) = 0, (9) 

the signs ± in Eqs. (|8|) and (|9|) being mutually independent. Obvious sym- 
metries of the equations allow us to fix only one set of the signs, as solutions 
corresponding to other combinations of the signs are actually the same. 

Getting back for a moment to Eq. (1) with ( = +1, it is worthy to 
note that, in this particular choice a = 1 and A = 1/V2 gives rise 

to an exact stationary solution interpolating between the real and imaginary 
asymptotic states, viz., 



(10) 



v(x) = (l/V2) [l + tanh (x/v 7 ^)] , w(x) = (1/V2) [l - tanh (x 
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However, if £ = — 1, no similar exact solution can be found. There is only 
one analytically tractable (and rather unphysical) case, which is A 1, i.e., 
a very broad ramp. In this case, defining £ = Ax, one can readily see that the 
^-derivative terms in Eqs. (0) and (0) are negligible, and an asymptotically 
exact stationary solution corresponding to A — > is 

w = 0, v 2 = atanh£ — 1 at tanh£ > 1/a, 

v = 0, w 2 = — atanh£ — 1 at tanh£ < —1/a, 

v = w = at |tanh£| < 1/a. (11) 

Of course, at small but finite A the fields which exactly vanish in various 
regions in the limit A = have small nonzero values at all £. 

In order to understand the structure of solutions to Eqs. (|) and (|7|), we 
consider the following auxiliary system, 



v X x — (« 2 v 2 + a 2 w 2 + 1) v + a g(x)v = 0, 
w xx — (a 2 v 2 + a 2 w 2 + 1) w — a g{x)w = 



(12) 



subject to b.c. 



v(—oo) = 0, v(oo) = 1, 
w(—oo) = 1, w(oo) = 0. 



(13) 



Here a = \J a — 1, and 




x e [rj, oo) 

g(x) = { x/r), x e [-77,77] , (14) 
x G (—00, —rj] 

with < 77 <C 1. The system ( |T2"D reduces to Eqs. if one substitutes 

g(x) by tanh(Ax), simultaneously replacing v and w by v/a and w/a in Eqs. 

It is easy to see that, due to symmetry arguments, v(x) = w(—x) and 
v(—x) = w(x), which in particular means that v(0) = w(0). Far from the 
origin, obvious outer-zone solutions to Eqs. (|T2| ) satisfying b.c. (0) are 

J v — il, w — 0, at x>rj, , v 

1 ii = 0, w = ±1, at x < —r/. 
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Inside the interval x G [—77, 77] one has, to the leading order in the small 
parameter 77, v xx = w xx = 0. An inner-zone solution to these equations, 
satisfying b.c. v{rf) = 1, v(—rj) = 0, 10(77) = and w(—rj) = 1, which are 
necessary for the continuous matching to the outer-zone solution (|l"5|)with 
the upper sign, is 

x 1 x 1 . . 

v = h - , w = h - • (16) 

2 77 2 ' 2 77 2 v ; 

Thus, to the leading order, a solution to the auxiliary equations (|I^)-([T3|), 
that we will call a type-A solution, is 

x G (-00, -rj\, 
v(x) = i (1 + x/t])/2, xG[— 77,77] 




x G [77,00), 

- x/t})/2, x G [-77,77], 
w(x) = \ 1, x G (—00, —77], 

x G [77, 00). 



(17) 



An example of the type-A solution is displayed in Fig. ([!]). 

We note that there also exists another inner-zone solution to the above 
leading-order equations v xx = w xx = 0, that satisfies b.c. v(rj) = — 1, 
v(— 77) = 0, w(rj) = and w(—rj) = —1: 

x 1 x 1 . . 

u = , w = . (18) 

2 77 2 ' 2 77 2 v ; 

This inner solution is matched to the outer one ([T^) with the lower sign. In 
order to make the latter solution also consistent with b.c. (|T3|) , we need to 
connect the two constant solutions v — ±1 in the outer region x > 77. To 
this end, we consider a region at x > rj where w = and g(a;) = 1. Then the 
second equation in the system fll2D is automatically satisfied, and the first 
equation yields 

v xx -(a-l) (v 3 -v) =0. (19) 

We choose an arbitrary point Xq > 2n and look for a solution to Eq. 
( [19]) in an interval (x — i],x + 77). If (a — 1) r/ 2 1, then the solution 
is 7; = (x — xo)/rj. If (a — 1) 7/ 2 = 0(1), then the solution to fll9D is v = 
tanh(w (a — l)/2 (x — xq)), whose expansion at small \x — xq\ is again (x — 
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xo)/rj. Thus, to the leading order, an approximate solution to Eqs. (|12|)-(|T3|) 
that we call a type-B solution, and which is valid for an arbitrary xq > 2rj, is 



v(x) 



o, 


X 


e 


(-00,-77], 


-(l+x/r])/2 


X 


e 




-1, 


X 


e 


[77, x - 77], 


(a; - x q )/t}, 


X 


e 


[x -r},x + r]], 


1 


X 


e 


([Xo + Tj, OO). 


1, 


X 


e 


(-00, -x - r/], 


-(x + x )/r], 


X 


e 


[-x - 77, -x 4 


-1, 


X 


e 


[— XO + 77, -77] 


(l-x/r,)/2, 


X 


e 




o, 


X 


e 


[77,00). 



(20) 



w(x) 



An example of this solution is displayed in Fig |2|. 

It is easy to see that the type-B solution (|20| ) can also be represented in 
the form 

v(x) = 2 f(-x - x ) - /(-x), w(x) = 2 f(x-x ) - f(x), (21) 
where /(x) is the f-component of the type-A solution, i.e., 



The representation 



(22) 



x e (-00, -77], 

x G [—77, 7]} 

x e [77,00), 

suggests that there exists a continuous family of 



0, 

(l + x/r/)/2, 
1, 



solutions to the auxiliary system ( 12]) , parametrized by xo, so that when 
xo = 0, the solution is of the type A, and when xo > 2e, the solution is 
of the type B. In the intermediate case xo G (0,277), the solution cannot be 
represented in the simple approximate form (plD , and can only be constructed 
numerically. In other words, the extra kinks in the v- and w- components, 
which distinguish the solutions of the type B from those of the type A, can 
be inserted, essentially, at any point, i.e., there is no unique "equilibrium" 
point at which the extra kinks have to be placed in order to provide for a 
stationary solution. 

Once the general structure of the solutions to the auxiliary system (|Tj) 
is understood, it is natural to check whether stationary solutions of the un- 
derlying equations (8) and (9) follow the same pattern. First, in order to 
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check the existence of type-A stationary solutions to Eqs. (8) and (9), we 
have performed numerical integration of the boundary- value problems || by 
means of the finite-difference method in a domain of the size Ax = 100 (as it 
will be seen from results displayed below, this size is quite sufficient to con- 
sider the system as being practically infinitely large). The number of spatial 
points was 8192, so that the spatial stepsize was 0.0122. 

The numerical problem was thus reduced to a system of algebraic equa- 
tions of the order 2 • (8192 — 1), which were solved by means of the Newton's 
method with tolerance 10~ 4 . As the first step, we looked for type-A solu- 
tions, taking, as initial trials for v(x) and w(x), functions C/(C + e~ 7 x ) and 
C/(C + e 7 x ) with different values of C and 7, which seem close to the type-A 
solutions sought for. The calculations have been made over a large range of 
values of the parameters from Eq. (|5p, a = 2 — 1000 and A = 0.1 — 10000. 
In all type-A solution has been found indeed. 

The next step was to check the existence of the type-B solutions. To this 
end, we employed two different tests. First, we computed type-B stationary 
solutions by means of the algorithm described above but with initial trial 
functions given by the expression (21) , where f(x) was replaced by the pre- 
viously found type-A solution, for the same values of a and A, and various 
values of xq. 

Following this way, in most cases tested we could indeed obtain type-B 
solutions corresponding to, virtually, arbitrary values of xq, in accord with the 
results of the analysis performed above. However, the result was inconclusive, 
due to numerical problems, for the most intriguing cases in which Xq is such 
that the minima of the initial trial functions (|2l|) , with f(x) replaced by the 
previously found type-A solution, were not close enough to —1 (those values 
of xq are close to 0). In those cases the sequence of approximations produced 
by the Newton's method did not converge to a solution sought for within the 
required tolerance in a reasonable amount of steps. 

The second method, which, simultaneously, gave a direct information 
about the dynamical stability of the stationary solutions, was based on sim- 
ulating the full system of PDEs (8) and (9), this time using the expression 
fl2ip , with f(x) replaced by the corresponding type-A solution found at the 
previous step, as an iinitial configuration. In doing this, we again took dif- 
ferent values of the parameter x i n Eq. (f2~T|). 

The numerical simulations of the full nonlinear PDEs were performed 
using a spectral method with the FFT transform based on 8192 spatial points, 
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and the modified Euler's method for the time evolution with the step 5t = 
3 • 10 -5 , in a time interval which was checked to be sufficient to obtain a 
stationary solution (the largest interval was < t < 30). 

The values of Xq were chosen, in all the corresponding to 400 A; 

spatial discretization points, k = 1, . . . , 6 (i.e., x = 4.883, 9.766, 14.648, 19.531, 24.414, 29.297). 
We have found that the solution, essentially, remains almost constant in time. 
In particular, the value of x at which v and w cross zero, while varying from 
v = — 1 to v = 1 and w = 1 to w = — 1 respectively, remains virtually con- 
stant. This numerical observations support the conjecture that the expression 
( |2~1~D with f(x) replaced by the corresponding type- A stationary solution is a 
fairly good approximation to the type-B solution of Eqs. @ and (0). 

This approximation has been found to further improve as both a and A 
in Eq. (]5|) increase. As an illustration, we show in Fig. [3] the graphs of the v- 
components of the Type-A and type-B solutions found for a = 2 and A = 10 
at the time moment t = 30. The B-type solution was obtained, in accord 
with what was said above, from the initial condition in the form of Eq. 
with f(x) replaced by the type-A solution that was found directly from the 
stationary ODEs. 



As well as in the case of the auxiliary system (|i~2|) and (H), the existence 
of the B-type solutions can be qualitatively understood. Indeed, it is clearly 
seen in Figs. [1| and |2| that the abrupt jump between positive and negative 
values of, say, v takes place when the w-component is practically equal to 
zero and, moreover, we may replace tanh (Ax) by +1 or —1, as the above- 
mentioned jump takes place far from the place where tanh (Ax) is essentially 
different from ±1. Thus, as a matter of fact, we are dealing with a very 
simple stationary GL equation, (|i~9"D, about which it is well known that it 
has a single stable nontrivial solution in the form of a kink, v = ±y/a — 1 
tanh (^/(o- l)/2ar), cf. Eq. (4). Therefore, the jump in the type-B solution 
exactly corresponds to this kink solution. Moreover, it is known too that 
more complicated stationary solutions to the real GL equation, containing 
several kinks and/or antikinks, do not exist, which explains why we have 
never obtained solutions with a larger number of jumps. 
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3 Stability of the stationary solutions 



Results presented in the previous section suggest that both the type-A and 
type-B solutions are dynamically stable as stationary solutions to the PDE 
system. To test the stability in a more systematic way, we simulated the non- 
stationary equations (H) and (|7|), adding to the (previously found) stationary 
solutions small initial perturbations in the form of an arbitrary combination 
of spatial harmonics compatible with the boundary conditions. 

Actually, we, first of all, solved the non-stationary equations (^]) and (0) 
which were linearized around the stationary solutions found in the previous 
section. In this approximation, all the small perturbations added to an ar- 
bitrary type-A stationary solution quickly decay. Thus, we conclude that all 
these solutions are unequivocally stable. 

The situation with the type-B solutions is more tricky: perturbations 
generally decay, but there remain some residual localized perturbation pulses, 
which show a clear trend to be stuck exactly at the spots where the jumps 
of the type-B solution are located, see Fig. f|. The fact that the residual 
perturbations of the v- and w-components accumulate at the point where, 
respectively, the v- and w-components of the unperturbed stationary solution 
make their jumps, is quite easy to understand, as at these points both the 
v and w fields in the unperturbed solution are nearly equal to zero, hence 
the nonlinear suppression of the perturbations does not locally take place, 
allowing the perturbations to survive. 

Further simulations, based on the full stationary equations (f|) and (0), 
rather than their linearized versions, demonstrate that, eventually, the resid- 
ual perturbation pulses seen in Fig. |](b) decay, but very slowly. This is 
natural too, as the above-mentioned mechanism disabling the nonlinear sup- 
pression of the perturbations is only an approximate one, and it cannot 
provide for their survival over infinitely long time. 

4 Conclusion 

In this work, we have introduced a parametrically driven gradient Ginzburg- 
Landau model, which is subcritical in the absence of the parametric drive. 
In the case when the parametric drive is uniform in space, a kink solution to 
this model is stable, on the contrary to the well-known instability of the same 
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solution in the GL model without the parametric drive. However, a really 
nontrivial situation takes places when the parametric drive is itself subject 
to the kink-like spatial modulation, so that it selects real and imaginary con- 
stant solutions at x — ±00. In this case, we have found stationary solutions 
numerically, and also analytically for a particular case. The solutions look 
as being of two different types, viz., a pair of kinks in the real and imaginary 
components, or the same with an extra kink inserted into each component. 
We have showed that, in fact, both solutions belong to the same continuous 
one-parameter family of solutions. Simulations show that the former solution 
is always stable, while the latter one admits a special type of small perturba- 
tions that remain virtually constant in time, rather than decaying or growing. 
Eventually, these special perturbations decay, but extremely slowly). 

We acknowledge valuable discussions with A. A. Nepomnyashchy and W. 
Zimmermann, and thank Pablo Funes for his help with the software im- 
plementation. One of the authors (B.A.M.) appreciates hospitality of the 
Max Planck Institute for Physics of Complex Systems (Dresden, Germany), 
where a part of this work was done. H.G.R. acknowledges support from the 
Fischbach Fellowship. 



11 



References 

[1] P. Coullet, J. Lega, and Y. Pomeau, Europhys. Lett. 15, 221 (1991); P. 
Coullet and K. Emilsson, Physica D 61, 119 (1991). 

[2] I.V. Barashenkov, Yu.S. Smirnov, and N.V. Alexeeva, Phys. Rev. E 57, 
2350 (1998); I.V. Barashenkov and E.V. Zemlyanaya, Phys. Rev. Lett. 
83, 2568 (1999). 

[3] R.E. Kelly and D. Pal, J. Fluid Mech. 86, 433 (1978); L.P Vozovoi 
and A.A. Nepomnyashchy, Prikl. Mat. Mekh. 43, 998 (1979) [English 
translation: J. Appl. Math. Mech. 43, 1080 (1979)]; G.V. Levina and A.A. 
Nepomnyashchy, Prikl. Mat. Mekh 47, 402 (1983) [English translation: 
J. Appl. Math. Mech. 47, 341 (1983)]; P. Coullet and P. Huerre, Physica 
D 23, 27 (1986). 

[4] D.A.S. Rees and D.S. Riley, J. Fluid Mech. 199, 133 (1989); R. Schmitz 
and W. Zimmermann, Phys. Rev. E 53, 5993 (1996). 

[5] G. Hartung, F.H. Busse, and I. Rehberg, Phys. Rev. Lett. 66, 2741 (1991). 

[6] B.A. Malomed and A.A. Nepomnyashchy, Europhys. Lett. 21, 195 (1993); 
Phys. Lett. 244, 92 (1998). 

[7] L. Kramer, E. Ben- Jacob, H. Brand and M.C. Cross, Phys. Rev. Lett. 
49, 1891 (1982); Y. Pomeau and S. Zaleski, J. Phys. (Paris) 44, L135 
(1983); H. Riecke, Europhys. Lett. 2, 1 (1986); R.B. Hoyle, Phys. Rev. E 
51, 310 (1995). 

[8] R.L. Burden and J.D. Faires. Numerical Analysis (PWS Publishing Com- 
pany: Boston, 1980). 



12 



Figure 1: A general scheme of the type- A solution. The corners of the kinks 
in the v and w components are rounded since the exact solution, unlike the 



approximation (17) , does not admit a discontinuity of the first derivatives. 
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Figure 2: A general scheme of the type-B solution. 
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Figure 3: The v- components of the type-A and type-B solutions to Eqs. (8) 
and (9). The latter solution was obtained using the initial condition in the 
form fl2lD with f(x) replaced by the type-A solution, that was found directly 
from the stationary version of Eqs. (8) and (9). The upper panel displays 
the solution as a whole, while the lower panel is a blowup of the region where 
the B-type solution makes its jump. In this figure, a = 2, A = 10, and the 
time moment shown is t — 30. 
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Figure 4: (a) A particular example of the initial perturbation added to the 
stationary solution of the type B (y- and w-components of this perturbation 
are equal), that gives rise to practically constant residual perturbation pulses, 
shown in the panel (b), which are pinned at the points where the type-B 
solution makes its jump. In this figure, a = 8.0, A = 2.0, and the panel (b) 
pertains to the time moment t — 15. 
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